Monte Carlo Minimization of the Higgs-Landau Potentials 



Jai Sam Kim/ J.C. Toledano,^ and P. Toledano'^ 
^Dept of Physics, Pohang University of Science and Technology, Pohang 790-784, S. Korea 

"^Ecole Poly technique, Rte de Saclay, F-91128, Palaiseau, France 
^ Grp de Physique theorique, Faculte des Sciences, 33 rue St-Leu, F-80000, Amiens, France 

(August 12, 1997) 

We propose a Monte Carlo algorithm to search for the boundary points of the orbit space which is 
important in determining the symmetry breaking directions in the Higgs potential and the Landau 
potential. Our algorithm is robust, efficient and generally applicable. We apply the method to the 
Landau potential of the d— wave abnormal superconductor. 
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I. INTRODUCTION 



There are many branches of physics where the symmetry breaking is important yet its analysis is highly non-triviaL 
The Higgs potential and the Coleman- Weinberg potential in the unification theories are difficult to minimize except 
for simple cases. Extremization of the superpotential in the supersymmetric unification theories is non-trivial. The 
Landau potential in many branches of condensed matter physics is often times tricky to handle. 

The Michel- Radicati conjecture |^ that a group invariant potential has extrema at a few critical strata corresponding 
to the maximal little groups was handy and utilized by many model builders (^,^. The group theoretical method 
allows one to locate the extrema from the list of maximal little groups without solving the coupled third degree 
extremum equations. It stirred up considerable research efforts. Slansky ||] made lengthy tables of little groups of 
many representations that could be used in the grand unification model building. Abud and Sartori |^ ellaborated 
the geometrical structure of the strata of an irreducible representation. Stokes and Hatch jsj made a complete set of 
tables for the isotropy groups of the irreducible representations of 230 space groups. However, Toledano's 1^ found 
counter-examples to the conjecture. Moreover, the conjecture is less useful for a reducible representation where strata 
of maximal little groups are oftentimes not one-dimensional critical strata. 

Some years ago one of the authors 1^,^ devised an efficient method for finding symmetry breaking directions in those 
problems. His method was practical and could be used for any representations. The key point in the method was the 
concept of the orbit parameters, which are defined to be the dimensionless ratios of (polynomial or non-polynomial) 
invariant functions to the unique isotropic quadratic invariant. For a reducible representation the ratio is taken over 
the product of quadratic invariants of each irreducible representation. Thus an orbit parameter is a kind of angular 
parameter. The space of such orbit parameters occupies a confined region, which is called the orbit space. Being 
defined out of invariants each point in the orbit space remains unchanged by the symmetry operations and has a 
definite residual symmetry associated with it, its little group, or sometimes called the isotropy group. Kim showed 
that the most general quartic Higgs or Landau potential has the absolute minimum at the most protruding point of 
the orbit space, which corresponds most likely to one of the critical strata in the case of an irreducible representation 
but not necessarily so in the case of a reducible representation. He showed that if the Michel-Radicati conjecture is 
to be valid the orbit space boundary has to be concave except at the cusps corresponding to the critical strata. Once 
the orbit space is constructed then the minimization of the potentials and the identification of symmetries are easy. 
One can easily identify the minimum point, simply by a look at the orbit space in the two dimensional case. 

The advantage of this method is that it finds not only the critical strata of maximal little groups but also non-critical 
strata of submaximal little groups and the phase diagram can be constructed easily. It is powerful if used with the 
extensive tables made by Slansky ||] and Hatch and Stokes It is most useful when the representation dimension 
is large and the orbit space dimension is two or three. 

It has been known [9UlO|] that the orbit space is a bounded object with a hierarchical structure. The points 
corresponding to the higher symmetries lie hierarchically closer to the surface of the orbit space and the interior 
points correspond to the unique lowest symmetry. The points of higher symmetries always form the boundaries of 
those of lower symmetries. The isolated points corresponding to the highest symmetries are zero dimensional objects 
and thus the derivatives of the orbit parameters with respect to the elements of the representation vector, sometimes 
called the carrier space vector, are zero. Thus the Michel-Radicati conjecture is a natural consequence of the orbit 
space structure and thus a group invariant function always has extrema on the critical strata [llll . There are additional 
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extrema which depend on the particular choice of the invariant potential. In order to determine the absolute minimum 
one has to compare all these local extrema. Thus one has to get the full orbit space boundary to locate it. 

Nobody showed that the orbit space boundary of an irreducible representation is concave except at the cusps. It 
is more likely that a high dimensional representation has many convex boundary portions as the number of available 
maximal little groups are limited. Moreover, the boundaries of an orbit space of a reducible representation have many 
convex boundaries . Nevertheless one can use this group theoretical property and systematically search for the little 
groups and the corresponding invariant directions to construct the orbit space ||^,|^ . Several examples were illustrated 
in 

One often times deals with the projected subspace of the full orbit space. With respect to the subspace the above 
hierarchy theorem does not necessarily hold. When the symmetry group is a non-simple group it is awkward to make 
a list of isotropy groups and the associated invariant directions. In this case the above theorem seems to be much 
less useful. Also when the dimension of the representation is large it is non-trivial to find the invariant directions 
corresponding to the little groups. 



II. MONTE CARLO SEARCH OF THE ORBIT SPACE 



Thus it is necessary to have an alternative method for building the orbit space. Simple Monte Carlo methods that 
were tried in the past were not efficient for obtaining the orbit space boundary. We have devised a more efficient 
Monte Carlo method for building the orbit space. It is applicable to any representation and robust. The algorithm 
uses the von Neumann rejection method Jl^ and the over-relaxation method. 

Suppose that we are searching for the boundary curves of a two-dimensional orbit space, {9, (p). Let us assume that 
the orbit parameters 9 and are functions of M variables, Xm- Then our algorithm is stated as follows: 

(1) Initially N points are sampled in a uniformly random way. That is, TV sets of Xm are sampled in a 
uniformly random way from the interval [-1, 1]. Then one computes the center of gravity, {9ct4>c)i 
of these N points, {9i, 

(2) For the N points, one computes the angle of the vector, {9i — 9c, (pi — (pc) and assigns an angular 
bin. Let the total number of bins be D. 

(3) One then identifies the outermost points along each of the D discrete directions. The 'outwardness' 
is defined to be the distance of the point from the center of gravity. 

(4) Starting from each of these selected points, a random walker with Af — legs makes moves by moving 
one of its legs, Xm- If the move increases the distance, it is accepted. If not, the random walker 
steps the leg in the opposite direction with a different stride. If the move increases the distance, it is 
accepted. Otherwise it is rejected and the old value of x,„ is retained. 

If the random walker jumps over into another angular bin, its position is compared with the current 
outermost point along that direction. The more distant of the two is selected as the new outermost 
point in that angular bin. 

(5) Step 4 is repeated for all Xm- 

(6) Step 4 and 5 are repeated along all the angular directions. 

(7) Step 4, 5 and 6 are repeated a finite times Rt or until no further appreciable improvements are made 
within the allowed number of trials. 

(8) The search of the orbit space boundary is completed. 

The circumference of the orbit space is discretized into D bins and the outermost point in each bin is sought after. 
In the random search for the boundary points the random walker of a selected bin is encouraged when it moves 
towards the surface. If a move x^ — *■ a^m + (5 is unsuccessful then another move Xm ^ Xm — uj ■ 6 is tried, where uj 
is an adjustable relaxation parameter. Even when the random walker strays into another direction its move is not 
discarded but compared with the current outermost point. In this way the random walker finds the outward direction 
in most attempts. Therefore the convergence rate is very fast. 
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III. APPLICATIONS 



Recently there have been some activities to search for anisotropic superconductivity It has been confirmed 
that a d-wave pairing actuaUy occurs [isl] . We have been investigating a model of an abnormal superconductor with 
the high temperature symmetry ||l6|, G = 0(3) x U{1) x T. For the spin-0 c?— wave gap functions the 3-dimensional 
symmetric traceless complex matrix, , is assumed for the order parameter. Under 0(3) the indices i,j transform 
like a three dimensional vector. The U{1) transforms $ — > e^^^ and <&* e~*'^$. The time reversal symmetry T 
transforms <i) ^ $* and vice versa. 

Then the most general Landau potential upto the fourth degree in $ can be written as a polynomial of the basic 
invariant polynomials. 

= alo + Pll + 71 Ti + 72r2 (1) 
Each basic invariant polynomial is defined as follows: 

^0 = 1>y1>*. (2) 

Ti = $,^$;fe T2^<^,,-^,kKin. (3) 

where repeated indices imply summation from 1 to 3. It is obvious that these are invariant under the symmetry group 
G. 

The orbit parameters, 9 and (j) are defined as: 

O^JI, '^^5 (4) 
Along the chosen direction {9, (p), the Landau potential has a directional minimum value, 

4 (/3 + 716' + 720) 

and the absolute minimum occurs along a direction where (/? + 71^ + 72?!') has the smallest value. 

By choosing an appropriate basis, we can always have one of the two real matrices in a diagonal form. So there are 
actually 7 independent components in '^u- We can define these 7 independent components like: 



^ki=\QbQ\+i\utw\ (6) 




We have run the program with the set of parameters: N — 2100, Rt = 7, D — 360, S — 0.125 and oj — 0.1. 
The acceptance rate was 0.460. The result of the run is as shown in Fig. 1. The orbit space, (0,0), is a triangle 
with vertices at Pl=(l. 0,0.0), P2=(0.5,0.5), and P3=(0.3,0.3). PI is associated with the set of variables, {a = l,b — 
— l,s = t = 0,u=l,v = w^O), and its little group is D^^ x e*" x T. P2 is overlapped by two cusps, P2 = (a = 
l,b=l,s — t^u = v = w~Q) corresponding to 0^(2) x T and P4 = {a = l,b = —1, s = t = u = v = w = 0) 
corresponding to D4^ii x e"^ x T. P3 is associated with {a — l,b = Q, s = Q,t — l,u — v — w — 0) where Q = —2 — V3 
corresponding to the little group D4^h x e*''/^ x T. 

The little groups of these points are maximal isotropy groups of the high temperature symmetry G. The spacial 
part of the isotropy groups of Pi, P3 and P4 are all but their cmbcddings in G are all different. The elements 
of these subgroups are listed, in terms of the symbols used in the International Tables for X-ray Crystallography, as 
follows: 

H{Pi) = {Ri, R2, Rq, Rio}, 

{i?3, i?4, -R12} X e*'^, 
{i?5. Re, i?i3, Rii} X T, 
{R7,Rs,Ri5,Rw} xe" xT 

H{P2)^0'{2)xT (8) 
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H{P3) = {Ri, R2, R5, Re, R9, Rio, Ri3,Ri4}, 

{i?3, i?4, i?7, i?8, i?ll, i?12, -R15, i?i6} X e^^^^ X T 



(9) 



H{P4) — {Ri, R2, R5, Rq, R9, Rio, Ri3, R14} ^ T, 

{R3, R4, Rj, Rg, Rii, R12, Ri5, Rig} x x T 

The Landau potential (1) truncated at the fourth degree can have no lower symmetry than those of these four 
points 10,1). The degeneracy at P2 can be lifted if we include the sixth degree terms, (see Fig. 2) 

It is interesting that the orbit space has a convex boundary portion in the subspace defined by (f) and ^ = 
(<i>jj$*j.<i>fci<i>;TO$^„$* -)//^. It is shown in Fig. 2. This implies that if we include the relevant sixth degree term in 
the Landau potential non-maximal symmetries may be observed at low temperatures depending on the stengths of 
the expansion coefficients. 



IV. DISCUSSIONS 



For a d— dimensional orbit space, the distance measure is defined in the same way, i.e., the distance from the 
random walker's position to the center of gravity. In three dimensions, if the whole range of the two angles, polar 
and azimuthal, are divided into 180 and 360 bins, then the number of bins alone is huge, 64800. For even higher 
dimensions binning the directions is out of the question. We will need an adaptive Monte Carlo method in this case. 
Namely, we allocate more random walkers in regions of interest. However, cusps can be found easily with only a few 
bins. 

There are several parameters that the user selects, the number of angular bins D, the initial number of points N, 
the number of repeats Rt, the random walker's stride 6, and the relaxation parameter lu. The optimum choice of 
these parameters can only be found empirically. 

We tested our method using several Higgs potentials and the Landau potential of '^He. First of all, the outcome 
depends on the choice of the center of gravity heavily. So we cannot sensibly give a performance report in a neat 
table of numbers. We found that N and Rt must be selected in proportion to the dimension M of the representation. 
A large value of D may yield sharper boundary curves but it need not be larger than 360. Cusps could be obtained 
very easily even with very small D. 

The random walker has to pass many singular curves on its march towards the surface. These curves are orbits of 
submaximal isotropy groups. They are on the boundary surface in the full orbit space and the Jacobian of the orbit 
parameters with respect to Xm vanishes along those curves 0. In the projected orbit space parts of the curves are 
buried inside. If (5 or w is too small the random walker may get stuck at the curves and may not be able to cross 
these singular curves. So a reasonably large number, 0.25 0.75, must be chosen as long as the rejection rate is not 
increased too much. 

Another point we want to emphasize is that the maximality conjecture is destined to be broken. In order to close 
an M— dimensional volume with concave boundaries we need at least (M + 1) cusps. For example, 3 cusps for a two 
dimensional volume, 4 for three, etc. If there are less than (M + 1) cusps then portions of the boundaries must be 
convex. For example, if we have two cusps for a two dimensional volume then we can make a banana shape at best 
instead of a concave triangle. For a Lie group the dimension of a representation is unlimited whereas the number 
of maximal isotropy groups is limited. Thus we have less than the minimal number of cusps and we will see convex 
boundaries, which means the breaking of the maximality conjecture. That is why we see a convex boundary in Fig. 
2. 

In conclusion, we have devised a general and robust Monte Carlo algorithm for finding the symmetry breaking 
directions. The method can be widely used both for the Higgs potential and the Landau potential. Applying it to 
the Landau theory of the d— wave superconductivity, we found that three four-fold symmetries are possible at low 
temperatures, which agrees with the experiment done on Tl2Ba2Cu06+5 pSj . 

We thank Prof. L. Michel for helpful discussions. This work was funded by Pohang University of Science & 
Technology and Ecole Polytechnique (France). 
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FIGURE CAPTIONS 

Fig. 1 A random sampling of the orbit space: {9, (f). 
Fig. 2 A random sampling of the orbit space: ((/),^). 
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FIG. 1. A random sampling of the orbit space: {0,4>)- 
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FIG. 2. A random sampling of the orbit space: (0,^). 
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